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SUMMARY 



Equations are developed for the design of double- 
panel, air--heated windshields capaole of preserving 
visihility during flight under icing conditions. 

The use of the equations is illustrated by calcula- 
tions for the design of a typical windshield. 



INTRODUCTION 

The investigation of the problems of ice prevention 
by the National Advisory Committee for Aeronautics has 
included a study of various means for preserving visibil- 
ity through the airplane windshield during icing condi- 
tions. Intensive development of a double-panel, air- 
heated windshield for the prevention of ice has been 
undertaken in thic connection at the Ames Aeronautical 
Laboratory of the National Advisory Committee for 
Aeronautics. Several windshields of this type have been 
constructed and tried repeatedly in icing conditions. As 
noted in reference 1« the device appears to be satisfac- 
tory and, since the completion of the tests, installations 
(which have given satisfactory service) have been made on 
both commercial and military a'ircraft. The present bulle- 
tin is concerned with the thermal and pressure relations 
which must be considered in the design of such a windshield. 



SYMBOLS 



The symbols used throughout this bulletin, their mean- 
ing, and the relation of units employed in their evaluation - • 
are as follows: 
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A cross- sect ional area of air gap normal to flow 
direction, sq ft 

b gap width, ft 

C thermal conductance of outer panel, Btu/(hr)(sq ft)(^F) 

Cvs specific heat of circulating air at constant pressure, 
^ Btu/(lt))(^F) 

D inside diameter of pipe, ft 

Dq equivalent duct diameter, ft 

d gap thickness, ft 

f friction factor, no units 

G weight" velocity, llD/(hr)(sq ft) 

g acceleration due to gravity, 41,8 X lo'^ , ft/(hr)(hr) 

h heat- transfer coefficient of surface, 
Btu/(hr)(sq ft)(^F) 

K thermal conductivity of wiiidshield material, 
Btu/(hr)(sq ft)(QF/in.) 

k therraal conductivity of fluid flowing., 
Btu/(hr)(sq ft)(^F/ft) 

I thickness of windshield material in direction of 
heat flow, in, 

m hydraulic radius, ft 

K gap length in direction of air flow, ft • 

P average absolute pressure of circulating air, 
■ Ib/sq ft 

q heat- transfer rate per unit area, Btu/(hr)(sq ft) 

E gas constant for air, 53 f t- lb/ ( lb) ( '^F abs- ) 

S surface area of windshield, sq ft 

T average temperature of circulating air, 
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Ti temperature of circulating air entering gap, 

tj^ inner- surf ace temperature of outer panel, .^1 

t^ out er- surf ace temperature of outer panel, 

Vg^^ average specific volume of circulating air, eu ft/lb 

W weight rate of circulating air, Ib/hr 

w weight rate of circulating air per foot of windshield 
width, lh/(hr) (ft ) 

|JL coefficient of viscosity of fluid .flowing, nD/(hr)(ft) 

AP pressure drop of circulating air through gap, Ib/sq ft 

AT temperature drop of circulating air through gap, *^F 

Av change in specific volume of air, cu ft/lh 

DEVELOPMENT OF DESIGN EQUATIONS 



Although continiied effort is being made to determine 
the convective heat- transfer coefficient from the surface 
of a windshield during icing conditions, reliable data 
have not yet been established. The need for the heat- 
transfer data for the outer surface is obviated by assum- 
ing a rate of heat flow through the panel and an outer- 
surface temperature which are based on the experience of 
extensive flight tests with several airplanes. The heat 
flow rate employed, as is suggested in reference 1, is 
1000 Btu per hour per square foot. In order to prevent 
ice formation, it is apparent that the out er- surf ace tem- 
perature must be above 32^ F and, therefore, a value of 
50° F is assumed which provides a small margin of safety. 

Practical design considerations limit the temperature 
of the circulating heated air which may' be employed. If 
the temperature is below 100° F, the pressure drop required 
to deliver 1000 Etu per hour per square foot to a surface, 
the outer temperature of which is 50° F> becomes greater 
than the available dynamic pressure of the air stream. 
While a blower may be employed to increase the available 
pressure drop, this complication may be avoided in most 
installations. The strength of glass, plastics, and other 
windshield materials is reduced at temperatures above 
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200^ F and, for this reason, the inlet temperature of the 
circulating air has "been limited to 200 J. The develop- 
ment of the pressure and temperature relations, therefore, 
has Deen based upon an assumed circulating-air temperature 
range of 100^ I to 200° 

A desirable simplification of the design problem can 
be achieved by the assumption that the thermal losses 
through the inner windshield panel to the cabin interior 
are small and may be neglected. Such an assumption has 
been made in the development of the. equations. Further- 
more, it also has been assumed that the dimensions of the 
windshield, the allowable pressure drop in the circulating 
air, and the physical properties of the windshield mate- 
rials are known and fixed. 

The heat-transfer rate per unit area through the 
windshield outer panel is given by equation (38), page 25, 
reference 2c In the symbols of the present bulletins this 
equation becomes 

q = C (ti - to) (1) 
Substituting the design assumptions, 

q = 1000 Btu/(hr)(sq ft) 
and to= 50° F, 

and transposing, the inner- surf ace temperature of the 
outer panel is given by 

1000 + 50 C 
^ ^ ~ 0 

Eqt'.ation (2) is plotted in figure 1 and shows the 
variation of t^ with C,. the outer-panel thermal 
conductance. 



If the outer panel is made of laminations of differ- 
ent materials, the thermal conductance for the panel is 
determined by 



in 



(3) 



where ''^/k^p ^g/^s* and ' so f..0T th are for successive 
laminations. It is recommended that -values of thermal 
conductivities (k) he obtained from the manufacturers of 
the windshield. materials. 

The temp.erature drop of the circulat ing-. air- in., 
passing through the gap is given hy 

■ ^ • . ' ;•• • AT = -i^- • i. ' 

For air at 150 F, the mean allowable temperature of. the 
circulating air (cp)" equals 0*24 Btu per pound per degree 

Fahrenheit and, assuming" q. equal'^ lOGO Btu per hour per 
square foot, the above equation becomes 

^ 1000 S ^ 4 167 'S . ..(4) 
0.24 W * w • • • • 

Equation (4) is plotted in figure 2 and shows AT 
as a function of 3i the' surf ace area of the windshield, 
for various values of W, the c ir culat ing- air weight 
rate necessary to supply the heat required for ice 
prevent ion. 

« 

The heat- t ransf er coefficient for the inner surface 
of the outer panel may be det ermined -by -.equat ion (ll), 
page 169, reference 2, developed for flow in pipes. 



^ = 0.0326 



For air in the temperature range of this analysis 
= 0.73; hence ■ ■ . • ; . : 

=0..03r£at" ' (6) 



But . . - ■ & ^ ?• 



Substituting in equation (.5) and transposing gives 



.6 



h = 0.02 



/ \0« 8 

D \\xkJ 



(6) 



Although equation (5) was developed for flow in pipes, 
it can "be applied, with reasonable accuracy, to small 
^•ectangular" duct s by the substitution of the equivalent . 
duct diameter (De) for the p ip e .diaciet er (D). By defini- 
tion, the equivalent duct diameter is equal to four timss 
the- hydraul Ic radius (m) and the hydraulic radius - §quai s 
the duct cross- sect ional area divided by the wetted per- 
imeter. Thus, if the gap thickness is d and the gap 
width is b, ■ . 



^ 4bd 

^ ^ 2b + 2d 



(7) 



Thickness: of gap will be small compared to width; there- 
fore only slight error -will .be involved in dropping the 
term »2d'» from the de^ominat or "'of "the fraction in equa^- 
tion (7). Then, 



and, 



ffl = 



'e _ 



2d 
2d d 



.4, ... 4.2 

Substituting this value of De in equation (5) gives 

.0.8 



(8) 



But 



Therefore, equation (8) "becomes: 

h ■= 0.01 ^ 



(9) 



and, evaluating the physical constants (k and p,) 
within the temperature range considered in this analysis, 
equation (9) hecomes: 



h = 



0,0032 w 



0«8 



(10) 
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Equation (lO) is plotted in figure 3 and shows the 
variation of h with ' d, the gap thickness, for various 
values of w, the circulating-air weight rate per unit 
gap width. 

In establishing the size of the gap and the quantity 
of the circulating air, consideration must he given to 
the allowable pressure drop through the gap, as well as 
to the thermal requirements. The pressure drop through 
the gap is calculated from equation (35), page 130, ref- 
erence 2o In the symbols of this bulletin, it becomes 

g 2gm 

The first fraction may be omitted for this analysis since 
Av, the change in specific volume of the air passing 
through the gap, is negligible. The pressure drop per 
foot of windshield length then becomes 



N- 2gm 



(11) 



The value of the friction factor (f) is given hy 
equation (4), page 111, reference 2, as 



f = 0.049 

'ifT 



Substituting this value of f in equation (ll) leads to 
the expression 

N Vd/ V 2gm J ^ ^ 

But G = - = — = ^ 

A bd d 

D = 2d 

and " n = ^ 

2 

1 . a 

Then, ^ = 5.56 X 10 ^- (12) 



8 



In equation (13), according to the universal gas 
laws the average specific volume of the circulating air 
is determined "by 



Assuming an average absolute pressure of S118 pounds per 
square foot and average absolut e t emperatur e of 150° plus 
460^ ]?, or 610° F, the above eauation becomes 



V = = 15,3 cubic feet per pound 

2118 ^ 

Substituting this value , of Vg^^ in equation (.13) leads 
to the expression 

4P 8.51 x 10 (14) 

K d^ 

Equation (14) is plotted in figure 4 and shows the 
variation of Ap/n with d for various values of w. 
It is of importance to note that figure 4 applies to sea- 
level conditions only, since the value of used to 
evaluate v^^ was 2118 pounds per square foot, absolute. 

In order to find Ap/n at anv other design altitude, 
multiply the value of AP/k obtained in figure 4 by the 
ratio 2118/P; where P equals the absolute static pres- 
sure at design altitude, pounds per square foot. 



APPLICATIOE OP DESIM EQUATIONS. 



Before the design equations and the relations shown 
by the curves of figures 1 to 4 can be applied to the 
solution of a windshield design problem, certain data 
relating to the materials and dimensions of the windshield, 
and to temperature, pressure, and quantity of the circulat- 
ing air, must be known or assumed, Por a typical case 
these data are as follows^ 

1. Materials of outer panel: Assumed to be lucite 

between two sheets of glass 

2. Windshield dimensions: Width (equals b, gap width) 

assumed to be 1.5 feet, and length (equals K, 



gap length) assumed to te 2.0 feet. Thickness 
of lucite (l^) assumed to be O0O3I inch, and 

thickness of each sheet of glass (Ig) assumed 

to be O0I25 inch 

3. Thermal conductivity of materials in outer panel; 

For lucite, is 1.25 Btu/(hr)(sq ft)(^F/in.) 

and for glass, Kg is 6,0 Btu/(hr)(5q ft)(^F/in, 

Allowable pressure drop through air gap (AP): 10 
pounds per square foot 

5. Available temperature of circulating air (T.): 

Assumed to be 200° F 

6. A tentative value of W, the weif:ht rate for the 

circulating air, must be assumede Since a 
v-alue of W is determined (through use of the. 
design equations and curves) which, for an 
acceptable solution, must agree with the ten- 
tative value, this may change as successive 
approximations of the final solution are made^ 

The problem t hen* r es olves itself into finding what 
value of d, in conjunction with the correct value of W 
(found by successive approximations) will provide a heat- 
transfer rate of 1000 Btu per hour per square foot through 
the windshield outer panel and maintain the outer surface 
at a temperature of 50° F without exceeding either the 
available circulating-air temperature (T^) or the allowable 

circulating air pressure drop (AP), The solution of the 
air-heated windshield design problem proceeds in the 
following manners 

lo Calculate C, using equation (3) and the known, 
or assumed, data regarding the thickness and 
thermal conductivity of the materialso 

2, Refer to figure 1 and, using the confputed value 

of 0, obtained t^^ 

3, From the known, or assumed, values of AP and 

determine AP/Ue 

Assuming a tentative circulating-air weight rate 
(Wi) determine w (equals W^^/b), a tentative 
value of the circulating-air weight rate per 
foot of windshield widthc 
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5. Refer to figure 4 and» using the computed values 

of Ap/h and w, obtain d. 

6. Kefer to figure 3 and* from* the values of d and 

w previously determined, obtain h. 

7. Prom h, t^, and q (equal to 1000 Btu per hour 

per square foot), calculate the average 
circulating-air temperature in the gap (Tgy) 

by the fundamental relation q equals h 
(Tav - *i)» or lav equals (1000 + hti)/h. 

8. With the assumption of a uniform temperature 

gradient in the gap, knowing Tg^^ and Tj_, 

calculate AT from the relation AT equals 
2(Ti - Tav)- 

9. Prom h and N, determine the surface area of 

the windshield (S). 

10. Refer to figure 2- and, from S and AT, obtain a 

second circulating-air weight rate (W^). 

11. Refer back to step 4 to see if the assumed value 

of agrees with the value obtained in step 

10, W-^ If the difference between W;^ and 
ffs is less than 10 percent of vV^ then the 
design may be considered complete; if not., 
then assume another weight rate (W3) halfway 
between b ■ and repeat procedure 

beginning with step 4. 

This process should be continued until the 
difference between the last assumed and final 
resultant values of W is less than 10 percent 
of the former. 

Por purposes of illustration, an example of an appli' 
cation of the de^sign procedure is presented, using the 
data assumed for a typical windshield installation. 



1 ^ __i ; 

Z, \Kg/ 1.25 \ 6 J 



= 15.1 Btu/(hr)(3q ft)(^P) 
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2. from figure 1, = 115° F 

3. Ap/n = 10/2 = 5 Ib/Csq ft) (ft) 

4. Assuming = 450 Ho/hr 

W, 450 , , . , 

w = -L = =- 300 l"b/Ur)(ft) 

"b 1.5 

5. From figure 4, d = 0.017 ft 

6. Prom figure 3, ii = 18.2 Btu/(hr)(sq ft)(°F) 

„ , 1000 + h ti 1000 + 18.2(115) ^ 0 

7. = i = -r.-. = 170 F 

h 18.2 

8. AT = 2(T. - T^) = 2(200 - 170) = 50° F 

9. S = "bN = 1 . 5 X 2 = 3 sq f t 

10. Prom fi§;ure 2, W = 210 Ib/hr 

11. Since the difference between Wj and is 

greater than 10 percent of , a second assump- 
tion will be made. As a second approximation, 
assume W3 to be halfway between 450 pounds 

per hour and 210 pounds per hour, or 330 pounds 
per hour and, returning to step 4, repeat the 
design process. 

4. w = ^^5. = 220 lb/(hr)(ft) 

1.5 

5. d = 0.014 ft 

6. h = 17 Btu/(hr)(sq ft)(°F) 

7. = 10Q0-±-l?(ll5) ^ ,740 J 

17 

8. AT = 2(200 ~ 174) = 52° F 



10. 



W4 = 2 50 Ib/hr 
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The difference between W^^ and W ^, 80 pounds per 
hour, is still too large. The process was continued and 
satisfactory agreement was reached at a \7eight rate of 
approximately 260 pounds per hour. The corresponding gap 
thickness was 0.013 foot (approximately 5/32 inch) and 
the tenrperature drop through the gap v?as 48^ F. 

After the windshield ' has heen installed, flight tests 
are usually made to determine the performance* Natural 
icing conditions are necessary for a final check, "but an 
estimate of the ice-'pr event ion capahilitj'^ of the installa- 
tion can he made from flights in non-icing conditions* 
If instrumentation can be installed to measure the outer- 
surface temperatures the adherence of the system to the 
design assumptions can he estahiiehed. If the surface 
temperature cannot he measured, a flight in 40° F a-ir will 
prohahly result in a windshield out er^ surf ace temperature 
above 50^ P, and a measured heat- transfer rate through the 
outer panel of more than 1000 Btu per hour per square foot 
indicates, under these conditiOiis> a satisfactory design. 
It should he rem.embered that the heat- transfer require- 
ment of 1000 Btu per hour per square foot of surface area 
is based \ipon the flight speeds of reference 1 and may. 
not be conservative above 20=0 mil^s per hoar. Further 
investigations of the design requirements are being car- 
ried out by AAl, using air-heated windshields under natural 
icing condit ions, . and should result in a better unde-rs-tand- 
ing of the problem* 
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Windshield surface area, S. sq f t 
Figure 2.- Circulating air temperature drop through the windshield air gap as 

a function of windshield surface area, at various circulating air 
weight rates, necessary to supply the heat required for ice-prevention. 
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Pigiure 3,- Variation of air-heated windshield outer-panel inner- 
surface heat-transfer coefficient with air gap thick- 
ness for various circulating air weight rates per unit gap width. 
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Circulating air presstire drop per foot of gap length, A P/N.lVCsq ft) (ft) 

Figure 4.- Variation of circulating^ air pressure drop per unit length of 

air path with air gap thic3cne8s in an air-heated windshield at 
sea level for various circulating air weight rates per \mit gap width. 



